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Abstract 

A search for the decay /i+ — > e+7, performed at PSI and based on data from 
the initial three months of operation of the MEG experiment, yields an upper 
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limit on the branching ratio of BR(^+ ^ e+7) < 2.8 x 10"" (90% C.L.). This 
corresponds to the measurement of positrons and photons from ~ 10^"* stopped 
/Lt "'"-decays by means of a superconducting positron spectrometer and a 900 htre 
hquid xenon photon detector. 

Keywords: Muon decay, lepton flavour violation. 
PACS: 13.35.Bv, 11.30.Hv 



1. Introduction 

We report here on the results of a search for the lepton flavour violating 
decay /x"*" — >■ e"'"7, based on data collected during the first three months period 
of the MEG experiment. This operates at the 590 MeV proton ring cyclotron 
facility of the Paul Scherrer Institut (PSI), in Switzerland. 

Lepton flavour conservation in the Standard Model (SM) is associated with 
neutrinos being massless. Recent observations of neutrino oscillations [l[ imply 
a non-zero mass and hence the mixing of lepton Savours. However, in mini- 
mal extensions to the SM, with finite but tiny masses, charged lepton flavour 
violating processes are strongly suppressed and beyond experimental reach. 

Additional sources of lepton flavour violation (LFV) [2|, [sl |j] appear in the- 
ories of supersymmetry, grand unification or in extra dimensions, giving pre- 
dictions that have now become accessible experimentally. Hence, the present 
lack of observation of a signature of charged LFV may change with improved 
searches and reveal new physics beyond the SM or significantly constrain the 
parameter space of such extensions. 

The strongest bounds on charged LFV come from the muon system, with the 
current limit for the branching ratio BR(^"^ e"^7) < 1.2 x 10"" (90% C.L.), 
set by the MEGA experiment 

2. Experimental Principle 

The — > e+7 process is characterized by a simple two-body final state, 
with the positron and photon being coincident in time and emitted back-to- 
back in the rest frame of the muon, each with an energy equal to half that of 
the muon mass. 

There are two major sources of background, one from radiative muon decay 
(RMD) ^+ — > e^VeV^'^ and the other from accidental coincidences between a 
high energy positron from the normal muon decay /x"*" e'^VeD^^ (Michel decay) 
and a high energy photon from sources such as RMD, positron annihilation- in- 
flight or bremsstrahlung. Both types of background can mimic a ^ ^ ej event 
by having an almost back-to-back photon and positron. It can be shown 0, 
taking into account the muon rate as well as the acceptances and resolutions, 
that the accidental case dominates. 

Hence the key to suppressing such backgrounds lies in having a continuous 
muon beam, a good quality beam transport system and precision detectors with 
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excellent spatial, temporal and energy resolutions. This is the basis for the novel 
design of MEG. 

3. Experimental Layout and the MEG Detector 

A schematic of the experiment is shown in Figure [T] Surface muons of 
28 MeV/c from one of the world's most intense sources, the ttES channel at PSI, 
are stopped in a thin, partially depolarizing polyethylene target, placed at the 
centre of the positron spectrometer. To facilitate a stopping rate of 3 x lO^/i+s"""^ 
in the 18 mg/cm^ thick target, with minimum beam-related background, a Wien 
filter and a superconducting transport solenoid (BTS) with a central degrader 
system are employed. The MEG beam transport system so, cleanly separates 
{7.5a) the eight times higher positron contamination to provide a pure muon 
beam. The use of a helium environment in the spectrometer, together with a 
slanted target, ensures minimal multiple scattering for both the muons and the 
out-going positrons. This is essential also for limiting background production, 
such as annihilation-in-flight and bremsstrahlung, in the acceptance region, cen- 
tred around 90° to the incoming beam. 

Positrons originating from muon decay are analyzed in the COBRA (COnstant- 
Bending-RAdius) spectrometer consisting of a thin- walled superconducting mag- 
net with a gradient magnetic field, a tracking system of low-mass drift chambers 
and two fast scintillator timing-counter arrays. 

The gradient magnetic field in the spectrometer, ranging from 1.27 Tesla at 
the centre to 0.49 Tesla at either end, is designed such that positrons emitted 
from the target with the same momentum follow trajectories with an almost 
constant projected bending radius, independent of their emission angle. This 
allows a preferential acceptance of higher momentum particles in the drift cham- 
bers as well as sweeping away particles more efficiently, compared to a uniform 
field. 

The drift chamber system (DCH) consists of 16 radially aligned modules, 
spaced at 10.5° intervals, forming a half-circle around the target. Each drift 
chamber module contains two staggered layers of anode wire planes each of nine 
drift cells. The two layers are separated and also enclosed by 12.5 /xm thick 
cathode foils with a Vernier pattern structure. The chambers are operated with 
a helium:ethane (50:50) gas mixture, allowing this low-mass construction to 
total 2.0 X 10^"^ Xq along the positron trajectory. 

Positron timing information originates from fast, scintillator timing-counter 
arrays (TO), placed at each end of the spectrometer. Each array consists of 15 
BC404 plastic scintillator bars, with 128 orthogonally placed BCF-20 scintillat- 
ing fibres. Each bar is read-out at either end by a fine-mesh photomultiplier 
tube, while the fibres are viewed by avalanche photo-diodes. The precise timing 
and charge signals provide both the impact point on the TCs and give directional 
information on the positron. 

The photon detector is a 900 litre homogeneous volume of liquid xenon (LXe) 
that subtends a solid-angle acceptance of ~ 10%. It uses scintillation light to 



3 




Figure 1: Schematic of the MEG experiment, a) details of the COBRA positron spectrome- 
ter and LXe detector, showing the positron tracking chambers (DC), the scintillator timing 
counter arrays (TC) and the superconducting gradient-field solenoid COBRA, b) The MEG 
beam transport system, surface muons enter from the left. Also shown are the crossed-field 
separator (Wien filter), the superconducting transport solenoid (BTS) and the COBRA mag- 
net with the central thin stopping target surrounded by the various detectors. 



measure the total energy released by the 7-ray as well as the position and time 
of its first interaction. In total, 846 photomultiplier tubes (PMTs), internally 
mounted on all surfaces and submerged in the xenon, are used. The advantage 
of using liquid xenon is its fast response, large light yield and short radiation 
length. Stringent control of contaminants is necessary since the vacuum ultra- 
violet (VUV) scintillation light is easily absorbed by water and oxygen even 
at sub-ppm levels. The xenon is therefore circulated in liquid phase through 
a series of purification cartridges, and in gas phase through a heated getter. 
Both the optical properties of xenon as well as the PMT gains and quantum 
efficiencies are constantly monitored by means of LEDs and pointlike ^"^^Am a- 
sources deposited on thin wires stretched inside the active volume. The detector 
is maintained at 165 K by means of a pulse-tube refrigerator with a cooling power 
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of 200 W. 

To select matched photon and positron candidates in a high rate, continuous 
beam environment and store sufficient information for offline analysis requires 
a well matched system of front-end electronics, trigger processors and data ac- 
quisition (DAQ) software. 

The front-end electronics signals (2748) are actively split and go to both 
the trigger and the in-housc designed waveform digitizer boards. The latter are 
based on the multi-GHz domino ring sampler chip (DRS), which can sample ten 
analogue input channels into 1024 sampling cells each at speeds of up to 4.5 GHz. 
The sampling speed for the drift chamber anode and cathode signals is 500 MHz, 
while that of the PMT signals from the photon detector and timing counters 
is 1.6 GHz. This strategy gives maximum flexibility, allowing various read-out 
schemes, such as zero suppression, on an cvent-by-cvent basis for various trig- 
ger types. The system achieves an excellent pile-up recognition, together with 
superior timing and amplitude resolutions, compared to conventional schemes. 

The trigger is based on fast information from the two detectors using PMTs: 
the liquid xenon photon detector and the positron timing counters. It makes 
use of a subset of the kinematic observables from /x-decay at rest, requiring an 
energy deposit in the photon detector in an interval around 52.8 McV, a time co- 
incident positron hit on the timing counters and a rough collinearity of the two 
particles, based on their hit topology. The decay kinematics is reconstructed 
by electronics boards arranged in a triple layer tree-structure. The signal dig- 
itization is executed by means of a 100 MHz, 10-bit flash analogue-to-digital 
converter. A pre-scaled, multi-trigger event scheme is used for data-taking al- 
lowing calibration, background and signal events to be read-out together. The 
typical signal event rate was 5 Hz, and the total DAQ rate was 6.5 Hz, with an 
average livetime of 84%. 

In total, nine front-end compiiters are used for the DAQ, each sending an 
event fragment to a central event building computer over a Gigabit Ethernet 
link. An integrated slow-control system enables both equipment control and 
monitoring. 

A detailed GEANT 3.21 based Monte Carlo simulation of the full apparatus 
(transport system and detector) was developed and used throughout the exper- 
iment, from the design and optimization of all sub-systems to the calculation of 
acceptances and efficiencies. 

4. Monitoring £tnd Calibrations 

The long term stability of the MEG experiment is an essential ingredient in 
obtaining high quality data over extended measurement periods. Continuous 
monitoring and frequent calibrations are a prerequisite. Apart from such items 
as the liquid xenon temperature and pressure, the drift chambers gas composi- 
tion and pressure and the electronics temperature, a number of additional mea- 
surements must be performed to keep the subdetectors calibrated and synchro- 
nized. The three most important are nuclear reactions from a Cockcroft- Walton 
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(CW) accelerator, radiative muon decay (RMD) runs and pion charge-exchange 
(CEX) reaction runs. 

Three times a week during normal data-taking, 7-rays of moderate energy 
coming from nuclear reactions of protons on a Li2B407 target are used. Protons 
of variable energy (400 < Tp < 1000 keV) are produced by a dedicated CW- 
accelerator placed downstream of the experiment. The muon stopping target 
is automatically replaced by a remotely extendable beam pipe which places 
the nuclear target at the centre of the detector. Photons of E-y = 17.67 MeV 
from '^'Li{p, 7)^Be allow the monitoring of the LXe detector energy scale, while 
coincident 7's from ^-'^B(p, 7)^^C {E^^ = 4.4, 11.6 McV) detected simultaneously 
by the timing counter and the xenon detector allow the determination of time 
offsets of TC bars. 

Once a week an entire day of RMD data-taking at reduced beam intensity 
was performed, with the trigger requirements relaxed to include non back-to- 
back positron-photon events in a wider energy range. 

Two CEX runs {tt~p — > n'^n — > 7771) were also conducted, one at the be- 
ginning and one at the end of the data-taking period. Pion capture at rest in 
a liquid hydrogen target produces photons with energy 54.9 < < 83.0 MeV. 
By detecting one of these photons with the LXe detector and the other at 180° 
by means of a set of Nal crystals, preceded by a lead/scintillator sandwich, two 
mono-energetic calibration lines at the extremes of the energy spectrum are ob- 
tained. These enable measurement of the energy scale and uniformity. Dalitz 
decays (tt^ — )■ 7e"'"e~) were also collected by using a photon-positron coincidence 
trigger, and used to study the detector time synchronization and resolution. 

The combined use of all these methods enables the investigation of possible 
systematic variations of the apparatus. 

5. Event selection and resolutions 

The data sample analyzed here was collected between September and De- 
cember 2008 and corresponds to ~ 9.5 x 10^^ muons stopping in the target. At 
the first stage of the data processing, a data reduction is performed by selecting 
events with conservative criteria that require the time of the photon detector 
signal to be close to that of a timing counter hit, and at least one track to be 
detected by the drift chamber system. This reduces the data size to 16 % of the 
recorded events. The pre-selected data are again processed and those events 
falling into a pre-defined window (blinding-box) , containing the signal region on 
the 7-ray energy and the time difference between the 7-ray and the positron, 
are "hidden", i.e. written to a separate data-stream, in order to prevent any 
bias in the analysis procedure. Only the events outside the blinding-box are 
used for optimizing the analysis parameters and for studying the background. 

During the course of the data-taking period the light yield of the photon 
detector was c;ontinuously increasing due; to the purification of the liquid xenon, 
which was performed in parallel. Furthermore, an increasing number of drift 
chambers suffered from frequent high-voltage trips resulting in a reduction of the 
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positron detection efficiency by a factor of three over the period. The increase of 
the xenon hght yield was carefully monitored with the various calibration tools 
and it is taken into account in the determination of the energy scale. The trigger 
thresholds were also accordingly adjusted to guarantee a uniform efficiency. For 
the drift chamber system we adopted a normalization scheme which depends 
only on the ratio of the signal positron reconstruction efficiency relative to that 
of the Michel positron, in order to be insensitive to absolute efficiencies. 

A candidate /i"*" — > e^7 event is characterized by the measurement of five 
kinematic parameters: positron energy {Eg), photon energy (-E-y), relative time 
between the positron and photon (te-y) and opening angles between the two 
particles (^£7 and (f>e'y)- 

5.1. Positron energy, E^, 

The positron track is reconstructed with the Kalman filter technique p] , in 
order to take into account the effect of the multiple scattering and energy loss 
in the detector materials in the variable magnetic field. 

The positron energy scale and resolution are evaluated by fitting the kine- 
matic edge of the measured Michel positron energy spectrum at 52.8 MeV as 
shown in Figure [5J The fit function is formed by folding the theoretical Michel 
spectrum form with the energy-dependent detector efficiency, and the response 
function for mono-energetic positrons. The latter is extracted from the Monte 
Carlo simulation of — >■ 67 decays, and is well described by a triple Gaussian 
function (a sum of a core and two tail components). 

The resolutions extracted from the data are 374 keV, 1.06 MeV and 2.00 MeV 
in sigma for the core component and the two tails, with corresponding fractions 
of 60%, 33% and 7%, respectively. The uncertainty on these numbers is dom- 
inated by systematic effects and was determined by varying both the event 
selection and fitting criteria. 

5.2. Photon energy, Ej 

The energy calibration and resolution of the 7-ray at the signal energy is 
extracted from the two CEX periods. Figure [3] shows an example of the energy 
spectrum measured with 54.9 MeV photons from the CEX reaction. A small 
correction is made to take into account the different background present in the 
LXe volume during the operation with the pion beam. 

The line shape is asymmetric with a low energy tail due to 7-rays converting 
in front of the LXe sensitive volume. A 3D mapping of the parameters is also 
made, since they depend to some extent on the position of the 7-ray conversion, 
mainly on the conversion depth inside the detector (w). As an example, the 
average resolution for deep events {w > 2 cm) is measured to be AE/E = 
(5.8 ± 0.35)% FWHM with a right tail of ctr = (2.0 ± 0.15)%, where the 
error quoted includes the variation over the acceptance. The energy scale is 
constantly monitored by looking at the reconstructed 17.67 MeV energy peak 
from CW protons on Li and confirmed by a fit of the photon energy spectrum 
to the expected spectra from the /i^ — )■ e'^h'ei^fi'y decay, positron annihilation- 
in-flight and 7— ray pile-up, folded with the line-shape determined during the 
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Figure 2: Measured Michel positron energy spectrum. A solid line siiows the fitted function 
as described in the text. 

Tr'^-experiment. The systematic uncertainty on the energy scale is estimated, by 
a comparison of these measurements, to be < 0.4 %. 

5.3. Relative time, t^j 

The positron time measured by the scintihation counters is corrected by 
the time-of- flight of the positron from the target to the TC, as measured by 
the track-length in the spectrometer. The photon time is determined by the 
waveforms of the LXe PMTs and corrected by the line-of-flight that starts from 
the positron vertex on the target and ends at the reconstructed conversion point 
in the LXe detector. 

In Figure m the relative time distribution between the positron and the pho- 
ton in a normal physics run is shown; the RMD peak (outside of the blinding- 
box) is clearly visible above the accidental background. The ie^-peak is fit- 
ted in the region of 40 < i?-,. < 45 MeV and, by taking into account a small 
-E^-dependence observed in the 7r°-runs, the timing resolution for the signal is 
estimated to be at^^ = (148 ± 17) ps. The relative time between the LXe detec- 
tor and the TC was monitored over the whole data-taking period by observing 
the RMD time peak in runs at normal intensity, and was found to be stable to 
within 20 ps. 

5.4. Relative angles, O^-y and (jje-y 

The positron direction and decay vertex position are determined by project- 
ing the positron back to the target. The 7-ray direction is defined by the line 
linking its reconstructed conversion point in the LXe detector with the vertex 
of the candidate companion positron. The resolution of the angles between the 
two particles is evaluated by combining the angular resolution and the vertex 
position resolution in the positron detector and the position resolution in the 
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Figure 3; Measured energy spectrum for 54.9 MeV photons from a CEX run. 

photon detector. The positron angular resolution is evaluated by exploiting 
tracks that make two turns in the spectrometer, where each turn is treated as 
an independent track. The 9- and 0-resolution^ are extracted separately from 
the difference of the two track segments at the point of closest approach to the 
beam-axis and are ag = 18 mrad, = 10 mrad. Due to this difference, 9e^ 
and (j)e'y are treated separately in the analysis. The vertex position resolutions 
are measured, using the two-turn technique, to be ~ 3.2 mm and ~ 4.5 mm in 
the vertical and horizontal directions on the target plane respectively. These 
values were confirmed independently by a method which reconstructs the edges 
of several holes placed in the target. 

The position of the photon conversion point is reconstructed by using the 
distribution of the light seen by the PMTs near the incident position. The per- 
formance of the position reconstruction is evaluated by a Monte Carlo simulation 
and it is validated in a dedicated CEX experiment by placing a lead collimator 
in front of the photon detector. The average position resolutions along the two 
orthogonal front-face sides of the LXe detector and the depth direction (w) are 
estimated to be ~ 5 mm and ~ 6 mm respectively. 

On combining the individual resolutions, the averaged opening-angle resolu- 
tions of 21 and 14 mrad for 9ef and 9^67 are obtained respectively. 

6. Data analysis 

The analysis algorithms are calibrated and optimized by means of a large 
data sample in the side-bands outside of the blinding-box. The background level 



^taking the 2-axis as the beam-axis, 6 is defined as the polar angle, while (p is the azimuthal 
angle 
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Figure 4: The relative time distribution tej showing the RMD peak obtained during physics 
runs, for 40 < Ej < 45 MeV. The actual resolution function used in the analysis takes into 
account the different gamma energy. 



in the signal region can also be studied with the event distribution in the side- 
bands since the primary source of background in this experiment is accidental. 
The blinding-box is opened after completing the optimization of the analysis 



algorithms and the background study. The number of fr 



e+7 events is 



determined by means of a maximum likelihood fit in the analysis window region 



defined as 46 MeV < E. 
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< 60 MeV, 50 MeV < < 56 MeV, 

< lOOmrad. 
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< 1 ns. 



Oe-yl < lOOmrad and {(pe-y 

An extended likelihood function C is constructed as 
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where Ngig, -/Vrmd and A^bg are the number oi ^ ^ 67, RMD and acciden- 
tal background (BG) events, respectively, while 5, R and B are their respective 
probability density functions (PDFs). A^obs(= 1189) is defined as the total num- 
ber of events observed in the analysis window and N = iVgig-l-iVRMD+ABG- The 
signal PDF S is the product of the statistically independent PDFs for the five 
observables (E^, E^, te-y, Oe-f and 4>e-i), each defined by their corresponding de- 
tector response function with the measured resolutions, as previously described. 
The RMD PDF R is the product of the PDF for ^07, which is the same as that 
for the signal and the PDF for the other correlated observables {E^, E^, 6eq 
and 4>e-t)- The latter is formed by folding the theoretical RMD spectrum [6| 
with the detector response functions. The BG PDF B is the product of the 
background spectra for the five observables, which are precisely measured in 
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the data sample in the side-bands outside the bhnding-box. The position de- 
pendence of the resolutions in the case of the 7— ray is taken into account in the 
PDFs, together with all their proper normalizations. The event distributions of 
the five observables for all events in the analysis window are shown in Figure [SJ 
together with the projections of the fitted likelihood function. 




Figure 5: Projected distributions for each observable, containing all events in the analysis 
window. A solid line shows the likelihood functions fitted to the data. 

The 90 % confidence intervals on iVgig and ^Vrmd sltb determined by the 
Feldman-Cousins approach Q. A contour of 90% C.L. on the (A'^gig, A^rmd)- 
plane is constructed by means of a toy Monte Carlo simulation. On each point on 
the contour, 90 % of the simulated experiments give a likelihood ratio (£/>Cmax) 
larger than that of the ratio calculated for the data. The limit for TVsjg is 
obtained from the projection of the contour on the TVgig-axis. The obtained 
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upper limit at 90% C.L. is A^gig < 14.7, where the systematic error is included. 
The largest contributions to the systematic error are from the uncertainty of the 
selection of photon pile-up events, the photon energy scale, the response function 
of the positron energy and the positron angular resolution. The confidence 
intervals are calculated by three independent likelihood fitting tools, each with 
different schemes and algorithms. The results are all consistent. The expected 
number of RMD events in the analysis window is calculated to be 40±8, obtained 
by scaling the number of events in the peak of the ie-y-distribution, obtained with 
lower energy cuts, using the probability ratio in the PDFs. This expectation is 
consistent with the best estimate in the likelihood fitting of {25'^\q). 

The upper limit on BR(yK"'" — >■ e~^j) is calculated by normalizing the upper 
limit on A^^ig to the number of Michel positrons counted simultaneously with 
the signal and using the same analysis cuts, assuming BR(/i evv) w 1. 
This technique has the advantage of being independent of the instantaneous 
beam rate and is nearly insensitive to positron acceptance and efficiency factors 
associated with the DCH and TC detectors, as these differ only slightly between 
the signal and the normalization samples, due to small momentum dependent 
effects. The branching ratio can in fact be written as: 

1ST. fE trig .TC DCH i i 

where Nevp = 11414 is the number of detected Michel positrons with 50MeV < 
Ee < 56MeV; P — 10^ is the prcscale factor in the trigger used to select 
Michel positrons; /^^p = 0.101 ± 0.006 is the fraction of the Michel positron 
spectrum above 50 MeV; Cg^s/g^"! = 0.66 ±0.03 is the ratio of signal-to-Michel 
trigger efficiencies; A™/j4™ = 1.11 ± 0.02 is the ratio of signal-to-Michel 
DCH-TC matching efficiency; e^J^^/e^^-^ = 1.02 ± 0.005 is the ratio of signal- 
to-Michel DCEL reconstruction efficiency and acceptance; Af^ — 0.98 ± 0.005 is 
the geometrical acceptance for signal photons given an accepted signal positron; 
Ce-y = 0.63 ± 0.04 is the efficiency of photon reconstruction and selection crite- 
ria. The trigger efficiency ratio is different from one due to the imposition of 
stringent angle matching criteria at the trigger level. The main contributions to 
the photon inefficiency arc from conversions before the LXe active volume and 
selection criteria imposed to reject pile-up and cosmic ray events. 

The limit on the branching ratio of the /x+ — )• e+7 decay is therefore 

BR(^+ ^ e+7) < 2.8 x 10"" (90%C.L.) 

where the systematic uncertainty on the normalization is taken into account. 

The obtained upper limit can be compared with the branching ratio sensitiv- 
ity of the experiment with these data statistics. The sensitivity is defined as the 
upper limit of the branching ratio, averaged over an ensemble of experiments, 
which are simulated by means of a toy Monte Carlo, assuming a null signal and 
the same numbers of accidental background and RMD events as in the data. 
The branching ratio sensitivity in this case is estimated tobel.3xl0~^^, which 
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is comparable with the current branching ratio limit set by the MEGA exper- 
iment [5|. Given this branching ratio sensitivity, the probabihty to obtain the 
upper hmit greater than 2.8 x 10~^^ is 5 % if systematic uncertainties in the 
analysis are taken into account. 

7. Conclusion and Prospects 

A search for the lepton flavour violating decay /i^ e^7 was performed 
with a branching ratio sensitivity of 1.3 x 10"^"'^, using data taken during the 
first three months period of the MEG experiment in 2008. With this sensitivity, 
which is comparable with the current branching ratio limit set by the MEGA 
experiment, a blind likelihood analysis yields an upper limit on the branching 
ratio of BR(^+ ^ e+7) < 2.8 x 10"" (90% C.L.). 

The problem of the reduced performance of the drift chambers, due to high- 
voltagc trips, has been solved and the chambers functioned successfully during 
our 2009 run period. Additional maintenance to the LXe-detector has also 
resulted in a near optimal light-yield. Further improvements to the timing 
counter fibre detectors and the digitization electronics are in progress. 
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